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Abstract 

In the first part of this paper, we prove the existence of global strong solution for 
Korteweg system in one dimension. In the second part, motivated by the processes of 
vanishing capillarity-viscosity limit in order to select the physically relevant solutions 
for a hyperbolic system, we show that the global strong solution of the Korteweg 
system converges in the case of a 7 law for the pressure (P(p) = ap 1 , 7 > 1) to 
entropic solution of the compressible Euler equations. In particular it justifies that 
the Korteweg system is suitable for selecting the physical solutions in the case where 
the Euler system is strictly hyperbolic. The problem remains open for a Van der 
Waals pressure because in this case the system is not strictly hyperbolic and in 
particular the classical theory of Lax and Glimm (see [2UE]) can n °t be used. 

1 Introduction 

We are concerned with compressible fluids endowed with internal capillarity. The model 
we consider originates from the XlXth century work by Van der Waals and Korteweg 
[381122] and was actually derived in its modern form in the 1980s using the second gradient 
theory, see for instance [9"1 12CH 13?] . The first investigations begin with the Young-Laplace 
theory which claims that the phases are separated by a hypersurface and that the jump in 
the pressure across the hypersurface is proportional to the curvature of the hypersurface. 
The main difficulty consists in describing the location and the movement of the interfaces. 
Another major problem is to understand whether the interface behaves as a discontinuity 
in the state space (sharp interface) or whether the phase boundary corresponds to a 
more regular transition (diffuse interface, DI). The diffuse interface models have the 
advantage to consider only one set of equations in a single spatial domain (the density 
takes into account the different phases) which considerably simplifies the mathematical 
and numerical study (indeed in the case of sharp interfaces, we have to treat a problem 
with free boundary). 
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Let us consider a fluid of density p > 0, velocity field u € R, we are now interested in the 
following compressible capillary fluid model, which can be derived from a Cahn-Hilliard 
like free energy (see the pioneering work by J.- E. Dunn and J. Serrin in [9] and also in 
[H [12"| 117]). The conservation of mass and of momentum write: 

i £ + a,( P v) = o, 

-(p^) + d x {p%u e f) - ed x (p e d x u £ ) + d x (a( P y) = e 2 d x K, 
where the Korteweg tensor reads as following: 

divK = d x {p e Ktf)d xxP * + ~(K(p e ) +p e K(p e ))\d x p e \ 2 ) - d x (K( P e )(d xP e ) 2 ). (1.2) 

k is the coefficient of capillarity and is a regular function of the form n(p) = e 2 p a with 

2 

a G R. In the sequel we shall assume that K(p) = —. The term d x K allows to describe 
the variation of density at the interfaces between two phases, generally a mixture liquid- 
vapor. P = ap 1 with 7 > 1 is a general 7 law pressure term, e corresponds to the 
controlling parameter on the amplitude of the viscosity and of the capillarity. When we 
set v e = u e + ed x (\\i p e ), we can write (jl.ip on the following form (we refer to [13] for the 
computations): 

( 4p £ + ^(pV)-e^ e = 0, 
? (1.3) 
-(pV) + d x {p\u"W)) - ed x (p e d x v e ) + d x (a( P y) = 0, 

We now consider the Cauchy problem of (|1 .3j) when the fluid is away from vacuum. 
Namely, we shall study (|1.3p with the following initial data: 

p e (0,x) =pl{x) >0, u e (0,x) =u e (x), (1.4) 

such that: 

lim (p%(x),u%(x)) = (p + ~,u + '~), with > 0. 

We would like to study in the sequel the limit process of system (|1.3p when e goes to 
and to prove in particular that we obtain entropic solutions of the Euler system: 

d 

—p + d x (pv) = 0, 

q (1-5) 

i (pv) + d x (pv 2 ) + d x {apl) = 0, 

Let us now explain the interest of the capillary solutions for the hyperbolic systems of 
conservation laws. 

1.1 Viscosity capillarity processes of selection for the Euler system 

In addition of modeling a liquid-vapour mixture, the Korteweg also shows purely the- 
oretical interests consisting in the selection of the physically relevant solutions of the 
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Euler model (in particular when the system is not strictly hyperbolic). The typical case 
corresponds to a Van der Waals pressure: indeed in this case the system is not strictly 
hyperbolic in the elliptic region (which corresponds to the region where the phase change 
occurs) . 

In the adiabatic pressure framework (P(p) = p" 1 with 7 > 1), the system is strictly hy- 
perbolic and the theory is classical. More precisely we are able to solve the Riemann 
problem when the initial Heaviside data is small in the BV space. Indeed we are in the 
context of the well known Lax result as the system is also genuinely nonlinear (we refer 
to [21] ). It means we have existence of global C 1 -piecewise solutions which are unique in 
the class of the entropic solutions. 

This result as been extent by Glimm in the context of small initial data in the BV-space 
by using a numerical scheme and approximating the initial BV data by a (^-piecewise 
function (which implies to locally solve the Riemann problem via the Lax result). For 
the uniqueness of the solution we refer to the work of Bianchini and Bressan ([2]) who 
use a viscosity method. 

In the setting of the Van der Waals pressure, the existence of global solutions and the 
nature of physical relevant solutions remain completely open. Indeed the system is not 
strictly hyperbolic anymore. 

If we rewrite the compressible Euler system in Lagrangian coordinates by using the 
specific volume r = 1/p in (t, 00) and the velocity u, the system satisfies in (0, +00) x R 
the equations: 

d t r - d x u = 0, 
d t u - d x (P(r)) = 0, 

with the function P : 00) — > (0, 00) given by: 

P(r) = P(~), TG(ioo). 
The two eigenvalues of the system are: 



(1.6) 



\ 1 ( T ,v) = -yJ-P'(T), X 2 (t,v) = -^-P'(t). (1.7) 
The corresponding eigenvectors r%, r 2 are: 

Mr,v) = ( ) , Mt,v) = ( ) (1-8) 

Furthermore by calculus we obtain: 

P"(t) -P"(t) 
VX 1 (t,v)-w 1 (t,v) = — y , V\ 2 {t iV ) -w 2 {t,v) = — v ; (1.9) 



2^-P\r) 2\J-P'(t) 

We now recall the definition of a standard conservation law in the sense of Lax (it means 
entropy solutions): 

• The system is strictly entropic if the eigenvalues are distinct and real. 
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• The characteristics fields are genuinely nonlinear if we have for all (t, v), 

VAi (t, v) ■ W\ (t, v) 7^ and VA2 (r, v) • W2 (t, v) 7^ 0, 

for more details we refer to [33] . The definition of genuine nonlinearity is some kind 
of extension of the notion of convexity to vector-valued functions (in particular when 
we consider the specific case of the traveling waves). The previous assumptions aim at 
ensuring the existence and the uniqueness of the Riemann problem ( see |10j and |33j). 

When P is a Van der Waals pressure, we observe that the first conservation law 
([331 ED]) is far from being a standard hyperbolic system, indeed: 

• It is not hyperbolic (but elliptic) in (^-, X R, 

• the characteristic fields are not genuinely nonlinear in the hyperbolic part of the 
state space. 

Here the classical Lax-Glimm theory cannot be applied. In particular there doesn't exist 
any entropy-flux pair, which suggests that the entropy framework is not adapted for 
selecting the physically relevant solutions. In order to deal with this problem, Van der 
Waals and Korteweg began by considering the stationary problem with null velocity, and 
solving VP{p) = 0. For more details we refer to 3JJ. It consists in minimizing in the 
following admissible set 

A = {p € L\VL)/W{p) € L 1 ^), / p(x)dx = m}, 

Jq 

the following functionnal 

F[p] = [ W(p(x))dx. 
Jo. 

Unfortunately this minimization problem has an infinity of solutions, and many of them 
are physically irrelevant. In order to overcome this difficulty, Van der Waals in the XIX- 
th century was the first to regularize the previous functional by adding a quadratic term 
in the density gradient. More precisely he considered the following functional: 

Flocal= [ {W(p^x)) +1 j\Vf\ 2 )dx, 

with: 

Aocai = H\fi)nA . 

This variational problem has a unique solution and its limit (as e goes to zero) converge 
to a physical solution of the equilibrium problem for the Euler system with Van der Waals 
pressure, that was proved by Modica in [28J with the use of gamma-convergence. 

By the Euler-Lagrange principle, the minimization of the Van der Waals functional 
consists in solving the following stationary problem: 

VP(//) = 7 eV VA//, 

where the right-hand side can be expressed as the divergence of the capillarity tensor. 

Heuristically, we also hope that the process of vanishing capillarity-viscosity limit 
selects the physical relevant solutions as it does for the stationary system. This problem 
actually remains open. 



4 



1.2 Existence of global entropic solutions for Euler system 

Before presenting the results of this paper let us recall the results on this topic in these 
last decades. We shall focus on the case of a 7 pressure law P(p) = ap 1 with 7 > 1 
and a positive. Let us mention that these cases are the only ones well-known (essentially 
because the system is strictly hyperbolic in this case and that we can exhibit many 
entropy- flux pairs). Here the Lax-Glimm theory can be applied, however at the end of 
the 70's, one was interested in relaxing the conditions on the initial data by only assuming 
po and uq in L°° . 

In the beginning of the 80's Di Perna initiated this program, consisting in obtaining 
global entropic solutions for L°° initial data. 

Indeed in [8j , Di Perna prove the existence of global weak entropy solution of (|1.5p 
for 7 = 1 + 2^T[ and 7 = 2k + ^ + 1 (with k > 1), d > 2 by using the so-called 
"compensated compactness" introduced by Tartar in [35J. This result was extended by 
Chen in [I] in the case 7 G (1, |] and by Lions et al in [26] in the case 7 G [3, 00). In [25] . 
Lions et al generalize this result to the general case 7 G (1,3), and finally the case 7 = 1 
is treated by [IS]. We would like to mention that these results are obtained through a 
vanishing artificial viscosity on both density and velocity. 

The problem of vanishing physical viscosity limit of compressible Navier-Stokes equations 
to compressible Euler equations was until recently an open problem. However Chen and 
Perepelista in [5] proved that the solutions of the compressible Navier-Stokes system 
with constant viscosity coefficients converge to a entropic solution of the Euler system 
with finite energy. This result was extended in |19] to the case of viscosity coefficients 
depending on the density. 

Inspired by [5] and [19] , we would like to show that the solution of the Korteweg system 
(jl.3p converges to a entropic solution of the Euler system with finite energy when the 
pressure is a 7 law. To do this, we will prove for the first time up our knowledge the 
existence of global strong solution for the Korteweg system in one dimension in the case 
of Saint- Venant viscosity coefficients. By contrast, the problem of global strong solutions 
for compressible Navier-Stokes equations remains open (indeed one of the main difficulties 
consists in controlling the vacuum). This result justifies that the Korteweg system allows 
us to select the relevant physical solutions of the compressible Euler system at least when 
the pressure is adiabatic (P(p) = ap 1 with 7 > 1). The problem remains open in the 
case of a Van der Waals pressure. 

1.3 Results 

Let us now describe our main result. In the first theorem we prove the existence of global 
strong solution for the Korteweg system (|1.3|) . 

Theorem 1.1 Let p > 0. Assume that the initial data po and uq satisfy: 

< m < po < M < +oo,p -p£H 1 (R),v G H 1 (R) n L°° (R) . (1.10) 

Then there exists a global strong solution (p, v) of \1.3\) on R + x R such that for every 
T > 0: 

p-pe L°°(0,T, H x (R)),p€ L°°(0,T, L°°(M)), 

v G L°°(0, T, iT^R)) n L 2 (0, T, H 2 (R)) and v G L°°(0,T, L°°(M)). 
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Finally this solution is unique in the class of weak solutions satisfying the usual energy 
inequality. 

Remark 1 We would like to point out that the problem remains open in the case of the 
Saint- Venant system, which corresponds to system M.3\) without capillarity. 

In the following theorem, we are interested in proving the convergence of the global 
solutions of system to entropic solutions of the Euler system (|1.5p . 

Theorem 1.2 Let 7 > I and (p t ,v t ) with m e = p t v e be the global solution of the Cauchy 
problem hl.3\l with initial data (pq,Vq) as in theorem Then, when e — > 0, there exists 

a subsequence of (p e ,m e ) that converge almost everywhere to a finite entropy solution 
(p,pv) to the Cauchy problem 111.5)) with initial data (po,PoVo). 

Remark 1 We would like to point out that Lions et al in [25] had obtained the existence 
of global entropic solution for 7 > 1 by a viscosity vanishing process, and the consid- 
ered regularizing system was exactly the Korteweg system modulo the introduction of the 
effective velocity. 

One important basis of our problem for theorem ll.2l is the following compactness theorem 
established in [5]. 

Theorem 1.3 (Chen-Perepelitsa Let ip G Cq(R) ; (77^ , (r) be a weak entropy pair 
generated by if). Assume that the sequences (p e (x,t),v € (x,t)) defined on R x R + with 
m e = p e v e , satisfies the following conditions: 

1. For any —00 < a < b < +00 and all t > 0, it holds that: 

[ [ {p e V +l dxdr <C(t,a,b), (1.11) 

JO Ja 

where C(t) > is independent of e. 

2. For any compact set K C R, it holds that 

ft {(p e y +e + p e \v e \ 3 )dxdT<C{t,K), (1.12) 
Jo Jk 

where C(t,K) > is independent of e. 

3. The sequence of entropy dissipation measures 

i/Vi"> e )t + ?V) mt )« are compact in H~l{R 2 + ). (1.13) 

Then there is a subsequence of (p e ,m e ) (still denoted (p e ,m e )) and a pair of mea- 
surable functions (p, m) such that: 

(p e ,m e ) ->• (p,m), a.ease^-0. (1-14) 

Remark 2 We would like to recall that the estimate t!.12\) was first derived by Lions et 
al in 126^ by relying the moment lemma introduced by Perthame in I30\j . 

The paper is arranged as follows. In section [2] we recall some important results on the 
notion of entropy enrtopy-flux pair for Euler system and on the kinetic formulation of 
Lions et al in [26J . In section [3j we show theorem 11.11 and in the last section 14.11 we prove 
theorem 11.21 
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2 Mathematical tools 

Definition 2.1 A pair of functions (rj(p,v), H(p,v)) or (rj(p,m),q(p,m)) for m = pv, 
is called an entropy-entropy flux pair of system \1.3\) . if the following holds: 

[ V (p,v)] t + [H(p,v)] x = 0, 

for any smooth solution of 111.5)) . Furthermore (r](p, v) is called a weak entropy if: 

77(0, u) = 0, for any fixed v. 

Definition 2.2 An entropy rj(p,m) is convex if the Hessian V 2 n(p,m) is nonnegative 
definite in the region under consideration. 

Such rj satisfy the wave equation: 

dttri = 6 2 p^d xx r,. 

From [26], we obtain an explicit representation of any weak entropy (77, q) under the 
following form: 

rf(p,m)= I x(p,s-v)ip(s)ds, 

Jr (2.15) 

H^(p, m)= (9s + (1 - 0)u)x{p, s - v)^(s)ds, 
Jr 

where the kernel x is defined as follows: 

x {o,v) = [p 2e - v 2 ] x + , X = 3 ~ 7 N > --, and 9 = 2^1 } 

A\H, ) IF j + , 2(7-1) 2 2 

and here: 

=t x for t > 0, 
=0 for t < 0, 



Proposition 2.1 (see fMj ) 

For instance, when ip(s) = \s 2 , the entropy pair is the mechanical energy and the asso- 
ciated flux: 

2 3 
m m t 

n*(p,m) = — + e(p), q*(p,m) = — ^ + e (p), (2.16) 
where e(p) = represents the gas internal energy in physics. 

In the sequel we will work far away of the vacuum that it why we shall introduce equi- 
librium states such that we avoid the vacuum. Let (p(x),v(x)) be a pair of smooth 
monotone functions satisfying (p(x),v(x)) = (p~' + ,v~' + ) when — h x > Lq for some 
large Lq > 0. The total mechanical energy for (|1.3|) in M with respect to the pair of 
reference function (p(x), v(x)) is: 

E[p,v](t)= [ (\p(t,x)\v(t,x)-v(x)\ 2 + e*(p(t,x),p(x)))dx (2.17) 
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where e*(p, p) = e(p) e (p) — e (p(p — p) > 0. The total mechanical energy for system (??) 
with k(p) = j is: 

E 1 [p,u](t) = [ (lp(t,x)\u(t,x)-u(x)\ 2 + e*(p(t,x),p(x)) + e 2 (d x p^) 2 )dx (2.18) 

JR z 

and the total mechanical energy for system (|1.3p is: 

E 2 [p,v](t)= I (lp(t,x)\v(t,x)-v(x)\ 2 + e*(p(t,x),p(x)))dx (2.19) 

JR z 

Definition 2.3 Let (po,vo) be given initial data with finite- energy with respect to the end 
states: (p^jV^) at infinity, and E[pq,vq] < Eq < +oo. A pair of measurable functions 
(p, u) : — > is called a finite-energy entropy solution of the Cauchy problem l\1.5\) 
if the following properties hold: 

1. The total energy is bounded in time such that there exists a bounded function 
C(E,t), defined on M + x R + and continuous in t for each E € M + with for a.e 
t > 0: 

E[p,v](t) <C(E ,t). 

2. The entropy inequality: 

^(p,v) t +q^(p,v) x <0, 
is satisfied in the sense of distributions for all test functions ip(s) 6 {±l,±s,s 2 }. 

3. The initial data {pq,vq) are obtained in the sense of distributions. 

We now give our main conditions on the initial data f j 1 . 4 j) . which is inspired from [S]. 

Definition 2.4 Let (p(x),v(x)) be some pair of smooth monotone functions satisfying 
(p(x),v(x)) = (p~ ,+ ,v~' + ) when — \-x > Lq for some large Lq > 0. For positive constant 
Co, C\ and C2 independent of e, we say that the initial data (Pq,Vq) satisfy the condition 
% if they verify the following properties: 

• Po>°> IrPo( x )\ u o( x ) ~ u( x )\ < Co < +00, 

• The energy is finite: 

I {\po( x M(x) ~ v(x)\ 2 + e*(p e (x), p(x)))dx <d< +00, 

JR z 

Jr Po( x ) 

In this section, we would like to recall some properties on the pair of entropy for the 
system (|l,5p . Smooth solutions of (|1.5p satisfy the conservation laws: 

d t r,(p,u) + d x H(p,u) = 0, 
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if and only if: 

Vpp = ^-y^-Vuu- (2-20) 

We supplement the equation 12.201 bv giving initial conditions: 

r/(0,u) = 0,77,(0, u) = #*). (2.21) 

We are now going to give a sequel of proposition on the properties of rj, we refer to |26] 
for more details. 

Proposition 2.2 For /)>0, n,w £R, 

• The fundamental solution of H2.20\) - [2.21\) is the solution corresponding to r/ p (0, u) = 
5(u) is given by: 

XM = (p^ l -ui 2 ) x + with A = ^=^T- (2-22) 

• The solution of \2.20\) - l2.21\) is given by: 

v (p, u ) = [ ^(Ox(p,£-«R, (2-23) 



• n is convex in (p, pu) for all p, u if and only if g is convex. 

• The entropy flux H associated with n is given by: 

H(p,u)= [ il>(S)[0t + (l-0)Z] X (p,£-u)dZ where 9 = ^. (2.24) 

JR z 

We now give a important result on the entropy pair (see [26] . lemma 4) . 

Proposition 2.3 Taking i^(s) = hs\s\, then there exists a positive constant C > 0, 
depending only on 7 > 1, such that the entropy pair {r ] i ',H^) satisfies: 

\^{p,u)\<(p\u\ 2 +pi), 

H^(p,u) > C-\p\u\ 3 + p^ +d ), for allp>0 andu£R, 

1 a (2.25) 

\ v i( P ,u)\<(p\u\+p e ), 

\vtm(p,u)\ < Cp' 1 . 

We are now going to give recent results on the entropy pair (r)^,q^) generated by tp G 
Cq(M) (we refer to [5] for more details). 

Proposition 2.4 For a C 2 function ip : R — > R, compactly supported on the interval 
[a, b], we have: 

supp(r]^), supp(q^) C {(p,m) = (p, pu) : u + p 9 > a,u — p 9 < b} : (2.26) 
Furthermore, there exists a constant Cw, such that, for any p > and «£K, we have: 



• For 7 G (1,3], 

\^(p,m)\ + \q^(p,m)\ < C^p. (2.27) 

• For 7 G (3, +oo), 

|r?^(p,m)| <C>, |^(p,m)|<C^(p + / +1 ). (2.28) 

• /f rf* is considered as a function of (p,m), m = pu then 

\r,i(p,m)\ + \pr,t m {p,m)\ < C*, (2.29) 
and, if rim is considered as a function of(p,u), then 

\vt(m,u)\ + Ip^nt^pu)] < (2.30) 

We now would like to express the kinetic formulation of (jl.5p introduced in ([26 J). 

Theorem 2.4 Let (p,pv) G L°°(R + , L 1 (R)) have finite energy and p > 0, then it is an 
entropy solution of 111.5]) if and only if there exists a non-positive bounded measure m on 
R + x R 2 such that the function x(Pi £ ~~ u ) satisfies: 

dtX + d x m+(l-9)u) X ]=d^m(t,x,0- (2.31) 



3 Proof of theorem 11.11 

We would like to start with recalling an important result due to Solonnikov (see [34|). 
Let po the initial density such that: 

< m < po < M < +00. (3.32) 
When the viscosity coefficient p(p) satisfies: 

M > c > for all/? > 0, (3.33) 
we have the existence of strong solution for small time. More exactly, we have: 

Proposition 3.5 Let (po,vq) satisfy \3. 32\) and assume that p satisfies \3. 33\) . then there 
exists Tq > depending on tuq, Mq, \\pq — and \\vqWh 1 such that \1.3\) has a unique 

solution (p,v) on (0, Tq) satisfying: 

p-peL^iH^R), a t pGL 2 ((0,T!) xM), 
v G L 2 {0,T ll H 2 {R)), d t v G L 2 ((0,Ti) x R) 

for all T\ < Tq. 

Remark 3 The main point in this theorem is that the time of existence Tq depends only 
of the norms of pq which gives us a low bounds on Tq of the system il.3\) . 
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In view of this proposition, we see that if we introduce a truncated viscosity coefficient 
Pn{p)- 

1 

p n (p) = max(/3, - ), 
n 

then there exists approximated solutions (p n ,v n ) defined for small time (0, To) of the 
system (|1.3p . In order to prove theorem 11.11 , W6 only have to show that (/?n> 

v n ) satisfies 

the following bounds uniformly with respect to n and T large: 

< m < Pn < M < +oo, Vt € [0,T], 

Pn-peL^iH^R)), (3.34) 



We are going to follow the method of Lions et al in [25J, indeed the main point is to 
prove that we can extend the notion of Riemann invariant or more precisely the kinetic 
formulation of proposition 12.41 to the system ()1.3[) . We recall that system (|1.3p has the 
following form: 



d 

T^Pn + d x (p n v n ) - ed xx p n = 0, 
d 

^ -^.(PnVn) + d x (p n v n v n ) - ed x (d x p n v n ) - ed x (p n d x v n ) + d x (a(p n y = 0, 

and we have finally: 
d 

777 Pn + d x (p n v n ) - ed xx p n = 0, 
at 

d 

-TriPnVn) + d x (p n v n v n ) - ed x d x (p n v n ) + d x (a(p n y = 0, 
at 



(3.35) 



(3.36) 



Following [25J and setting m n = p n v n we have for any pair of entropy flux (ry(p, u), H(p, u)) 
defined by (|2.23[) and ()2.24p where rj is a convex function of (p n ,mn)- We write rj = 
fj(p n ,rn n ): 

d t rj + d x H = efj p d xx p n + efj m d xx m n , 

= td xx 7] - e(fi pp (d x p n ) 2 + 2f] pm (d x p n )(d x m n ) + fj mm (d x m n ) 2 ) . 

Here we define p n such that: 

Pn = f]p P (d x p n ) 2 + 2fj pm (d x p n )(d x m n ) + f) mm (d x m n ) 2 

By proposition 12. 2| we can check that p n > 0. We obtain then that: 

dtv(Pn,v n ) + d x H((p n , v n ) - efjpd xx p n < inlx (0, +oo). 

By applying the same method than for proving the theorem 12 we obtain the following 
kinetic formulation: 

dtX + d x ([0£ + (1 - 6)v n ] X ) ~ d xxX = %m n on R 2 x (0, +oo), (3.37) 
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where fh n is a nonpositive bounded measure on M. 2 x (0, +00). Finally we recover the 
classical maximum principle by multiplying (|3.37p by the convex functions g(£) = (£ — 
£0)+ and g(£) = (£ — £0)- and integrating over M 2 x (0, +00). Indeed as we have that: 

— C < minf^o — p Q ) < max(t>o + p ) < C, 

X X 

and that: 

supp£ = [v - p ,v + p[0}. 
For ^0 large enough, we can show that: 

supp£o n suppx = 0. 

We have obtain then that: 

-C < mm(v - Pq) < v n - p 9 n < v n + p e n < max(u + p ) < C. 

X X 

In particular we obtained that p n and v n are uniformly bounded in L°°(0, T n , L°°(M)) or: 

sup (\p n (t,x)\ + \v n (t,x)\) < C , (3.38) 

xeM,tG(0,T n ) 



4 Proof of theorem 11.2 



4.1 Uniform estimates for the solutions of (II. 3|) 

First we assume that (p € ,v € ) is the global solutions of Korteweg's equations (|1.3p con- 
structed in theorem 11.11 and satisfying: 

p e (t, x) > c e (t), for some c e (t) > 0, (4.39) 

and 

lim (p e ,v e )(x,t) = (p ± ,u ± ). (4.40) 

x— >±oo 

Here we are working around a non constant state (p, v) with: 

lim (p,v)(x,t) = (p ± ,u ± ). 

x— ¥±00 

It is a simple extension of theorem 11.11 Our goal is now to check the properties (|l.lip . 
(|1.12p and (|1.13|) in order to use the theorem 11.31 of Chen and Perepelista (see [S]) in 
order to prove the theorem 11.21 

For simplicity, throughout this section, we denote (p, v) = (p e ,v e ) and C > denote the 
constant independent of e. 

We start with recalling the inequality energy for system (|1.3p , indeed by the introduction 
of the effective velocity we obtain new entropies (see [13J). 

Lemma 1 Suppose that Ei[pq,uq] < Eq < +00 for some Eq > independent of e. It 
holds that: 

sup £?i[/j,w](r) + e / / pu\dxdr < C '(t) , (4.41) 
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and: 

ft 



sup E 2 [p,v}(T)+e [ [ pvldxdr + e [ [ p 1 ' 2 p 2 x dxdr < C(t), (4.42) 
o<r<t Jo Jr Jo Jr 

where C(t) depends on E$, t, p, and u but not on e. 

Proof: It suffices to writes the energy inequalities for system (|1.3|) and from system 
(|1.5p . More exactly we have: 

— / (rj(p,m) — n(p, pu)dx + e f pu 2 x dx = q(p~ ,m~~) — q(p + ,m + ), 



dt _ 

with the entropy pair: 



2 ,: i 

m , . , . m 



r){p,m) = — + e(p), q(p,m) = ^ + me (p), 
with e(p) = p 1 ■ Since we have: 

e(p,p)>p(p e -p 9 ) 2 ,0 = ^, 

we can classically bootstrap on the left hand-side the term q(p~,m~) — q(p + ,m + ). 
Remark 4 Since vacuum could occur in our solution, the inequality 

/ pu 2 dxdr < C(t), 
o Jr 

in h4-4^\ ) * s much weaker than the corresponding one in [5]. That is why lemma [21 will be 
more tricky to obtain. 

The following higher order integr ability estimate is crucial in compactness argument. 

Lemma 2 If the conditions of lemma{l\ hold, then for any — oo < a < b < +oo and all 
t > 0, it holds that: 

[ [ p^ +1 dxdr <C(t,a,b), (4.43) 

Jo Ja 

where C(t) > depends on Eq, a, b, 7, t, p, u but not on e. 

Remark 5 The proof follows the same ideas than in the case of compressible Navier- 
Stokes equations when we wish to obtain a gain of integrability on the density. We refer 
to fEffl for more details. The proof is also inspired from Huang et al in IT9\j . 

Proof. Choose u G C£°(R) such that: 

< lj(x) < 1, u)(x) = 1 for x € [a, b], and suppw = (a — 1, b + 1). 

By the momentum equation of (|1.3p and by localizing, we have 

{P(p)uj) x = -(puvu) x + (P(p) + puv)uj x - (pv) t uj + e(pv x uj) x - epv x uj x . (4.44) 
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Integrating (|4.44p with respect to spatial variable over (— oo,x), we obtain: 

/x rx 
pvLody) t + [(puv + P(p))uj x - epv x co x . (4.45) 
-oo J —oo 

Multiplying (|4,45p by pcj, we have 

pP(p)oj 2 = — p 2 uvoj 2 + ep 2 v x uj 2 — (pu> / pvujdy)t 

J —oo 

/•X PX 

- (pu) x uj( / puudy) + pu [{puv + P(p))uj x - epv x u x ]dx, 



—oo 



(4.46) 

--ep 2 v x u 2 - (pu; I pvudy)t - (puu I pvudy) x 



+ puuj x \ pvudy + pui I {{puv + P(p))u x — epv x u x ]dx, 



We now integrate (|4.46p over (0,t) x M and we get: 

f [ ap~ t+1 uj 2 dxdT = e f [ p 2 
Jo Jm. Jo Jm. 

p px ft f px 

+ / (pqoj / PqVq uidy)dx + I I (puu x / pvujdy)dxdr (4.47) 
Jm J -oo Jo Jm J-oo 



ft f f px 

" 2 v x uj 2 — / (poj / pvudy)dx 



o 



+ / J (puj I [{puv + P(p))u x - epv x uj x ]dx)dxdT. 

Let 



A = {x : p(t, x) > p}, where p = 2 max( / o+, p— ), (4.48) 
then we have the following estimates by (|4.42p : 

\A[ < °^ - = d(t). (4.49) 

By (|4.48p . for any (t, x) there exists a point xo = xo(t,x) such that \x — xo\ < d(t) and 
p(t, xo) = p. Here we choose (3 = > 0, 

supp xesupp(w) ep /3 (t,j;) < epf + supp xesupp(uj)nA ep f3 (t,x), 

< 2epP + supp a . esupp ( £j)nA |e/(i,a;) - ep^(t,x )\, 

i-x +d{t) 

< 2e/ + supp a;6supp(£j)nj4 / \{3\\epP~ 1 (t,x)p x \dx, 

J XQ—d(t) 

b+l+2d(t) r (4.50) 



<2e/+ / |/3|p 2/3 ~ 1 da;+ / Jp^fidx, 

Ja-l-2d(t) Jm. 
rb+ 

< C(t) + / 



'a-l-2d(i) 
rb+l+2d(t) 

p^dx, 

-2d{t) 



< C(t). 
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Using (|4.50p . Young inequalities and Holder's inequalities, the first term on the right 
hand side of (|4.47j) is treated as follows: 

t 

p 2 v x uj 2 dxdT 



o 



< —e I j p 3 u) 4 dxdr + — e / / pv 2 dxdr, 

< C(t) + e / / P 3 uj 2 dxdT, (4.51) 

Jo JR 

< C(t) + C(i) / / P A ~Pu?dxdT, 

Jo Jr 

< C(t) + 5 [ [ pi +1 oj 2 dxdT, 

Jo Jr 

Here we have used the fact that 7 > |. By lemma [T] and the Holder inequality, we obtain 



X P 

pvudy\ < / \pv\dy, 

-00 Jsupp(tj) 

,2 



(4.52) 



<(/ prfy) 5 (/ frfdy)* <C{t). 

supp(w) J supp(w) 



Then: 



(4.53) 



(pa; / p-uwtiy) <ix | + j / (poo; / PoVo^dy) dx | 

J-00 JR J-oo 

+ I / / {puui x I pvojdy)dxdr\ < C(t). 
Jo Jr J -00 

Similarly, we have: 

ft p rx 

(pu (puv + P{p))u x dy)dxdT\ <C(t), (4.54) 



and 

r t r 

(pa; / pv x uj x dy)dxdT\ 



< e| / / [pu I p\v x \\uj x \dy)dxdT\, 

Jo Jr Jr 

< e| / ( / pa;cte)( / pv 2 x dy+ / pa;^y)dr|, 

Jo Jr Jr Jr 



(4.55) 



< C(t). 

Substituting (|4.5ip . (|4.53p - (j4.55f) into f|4.47j) and noticing the smallness of 5, we proved 
lemma [2J 

Lemma 3 Suppose that (po(x),vo(x) satisfy the conditions in the lemmas [7J Further- 
more there exists Mq > independent of e, such that 

Po(x)\v (x) - v(x)\dx < M < +00, (4.56) 
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then for any compact set K C R, if ZioWs that: 

[ [ (p~ t+d + p\v\ 3 )dxdT <C{t,K), (4.57) 
Jo Jk 

where C(t,K) is independent of e. 

Remark 6 In order to prove the inequality ft4.57\ ), we will use the same ingredients than 
in ]26§ where this inequality was obtained for the first time. 



Proof. We are now working with the function ip of proposition 12.31 If we consider rjti 
as a function depending of (p, v), we have for all p > and 



(\ni v (p,v)\<C, 

For this weak entropy pair (rp 3 ,H^), we observe that: 

^(p,0)=^(p,0) = 0, H^p,0) = e -p i9+1 [ \sf[l-s 2 ]i 

and: 

vt(p,0)=Pp° with /3 = [ \s\[l - s 2 ] x + ds. 



(4.58) 



By Taylor formula, we have: 

r]*(p,m) = f3p e m + r(p,m), (4.59) 

with: 

r(p,m) < Cpv 2 , (4.60) 
for some constant C > 0. Now we introduce a new entropy pair (rj, H) such that, 

rj(p, m) = rj^(p, m — pv~), H(p, m) = H^(p, m — pv~) + v~n^(p, m — pv~), 

with m = pv which satisfies: 

f ff(p, m) = (3p e+1 (v -v') + r(p, p(v - v~)), 
1 r(p,p(v — v~)) < Cp(v — v~) 2 . 

Integrating (|1.3|) i x rj p + (jl.3|) 9 x rj m over (0, t) x (—00,2;), we have: 

(rf(p,m) -r)(p ,m ))dy+ / q*(p, p(v - v~)) + v~rf(p, p(v - v~))dr 

= tq*(p ,0)+e/ rj m pv x dT-e / {rj mu pv 2 x + rj mp pp x v x )dydT. 
Jo Jo J-00 

By using (|4.58p . we obtain: 

pt px pt p 

|e / / Vmupvldydr] < Ce / pv 2 x dy dr < C '(t) , (4.63) 
Jo J-00 Jo Jk 
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/ / Vm P PPxVxdydT\ < Ce / p e 1 p\p x v x \dy dr <C(t), 

>o J-oo Jo Jw ^ 

<Ce [ f pvldydr + Ce f [ p 7 " 2 p 2 x dy dr < C{t). 



Substituting (|4.63p and (|4.64[) into ([4.62p . then integrating over K and using (|2.25j) . we 
obtain: 

/ / P e+1 + p\v -v~\ 3 dxdr 

JO J K 

<C{t) + C / \r]*(p,p(v-v-)\dxdT + Ce / p\v\\v x \dxdr (4.65) 
Jo Jk Jo Jk 

+ Ce \ \ p 1+9 \v x \dxdr + 2 sup | / (/ eta(p(y, r), (pv){y, r))dy)dx\ . 

Jo Jk re[0,t] JK J-oo 

Applying lemma [TJ we have: 

f [ \r)*(p,p(v-v-)\dxdT<C(t). (4.66) 
Jo Jk 

By Holder's inequality and (|4.50|) . we get: 

el [ p 1+e \v x \dxdT <Ce [ [ pv 2 x dxdT + C e f [ p 1+2e dxdT, 
Jo Jk Jo Jk Jo Jk 



< C(t) + C{t) I I p d dxdr, 

Jo Jk 

< C(t). 

We have now: 

e [ [ p\v\\v x \dxdr < —e I I pv 2 dxdr + -e J f pv 2 dxdr, 
Jo Jk % Jo Jk 2 J J K 

< C(t). 



(4.67) 



(4.68) 



Now we are going to deal with the last term on the right hand side of f|4.65j) . (|1.3p implies 
that: 

(pv - pv~) t + {pv 2 + P(p) - puu~) x = e(pv x ) x . (4.69) 

Integrating (|4.69p over [0,t] x (— oo,x) for x G K, we get: 

p X ft 
p(v - pv~)dy = I po(v - pv~)dy - \ {pv 2 + P(p) - puu~ - P(p~)) 

+ e / pv x dr. 
Jo 
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Furthermore: 

v((p(y, T ), (pv)(y,T))dy\ 

oo 

<| r {v(ppv)-Pp 0+1 (v-v))dy\ + \ [ X (3p e+1 (v-v))dy\ 

J —oo J —oo 

< I / (r(pp(v - v))dy\ + | f (3(p 9 - (p-) e )p(v - v))dy\ (4.71) 

J —CO J —oo 

+ f3(p-) e \ f p(v-v))dy\, 

J —oo 

<C(t) + (3(p-) e \ f p(v-v))dy\. 

J —oo 

By using (|4.56p . lemma [1] and [21 (|4,70p and (|4.7ip we conclude the proof of the lemma. 
4.2 Hf l(Rl) Compactness 

In this section we are going to take profit of the uniform estimates obtained in the previous 
section in order to prove the following lemma, which gives the ^~*(M+)-compactness of 
the Korteweg solution sequence (p e ,v e ) on a entropy- entropy flux pair. 

Lemma 4 Let if) € C 2 (R), rfi ,H^) be a weak entropy pair generated by ift. Then for the 
solutions (p e ,v e ) with m e = p e v e of Korteweg system |773J the following sequence: 

rf(p e , m e ) t + q^{p\m e ) x are compact in Hf^(M. 2 + ) (4.72) 



Proof: Now we are going to prove the lemma. A direct computation on (jl.3|) i x 

■qf(p e ,m e ) + (fL3])9 x rjti(p e ,m e ) gives: 

^(p%m e ) 4 +^^p%m e ) x . = 6(^(p%m e )(^K)-6r ? L(p%^)(/)(<.) 2 

Let K C M be compact, using proposition I2.4I (|2.3U|) and Holder inequality, we get: 
t r 



JO JK 

ft r rt 



(4.74) 



<Ce[ ! (p*)\(vl) 2 dxdT + Ce f [ ( P y- 2 (pl) 2 dxdr 
Jo Jk Jo Jk 

< C{t). 
This shows that: 

- ert u {p e M)p\vl) 2 - erit u {p e M)p'v%p% are bounded in L\[0,T] x K), (4.75) 
and thus it is compact in W,~ c ' pi (K+), for 1 < pi < 2. Moreover we observe that 

k/L(/,/>v)|<c^, 
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(4.76) 



then we obtain: 

JK 

f* f 4, ,4, ,4 

< / / €3 |p £ |3 3c?X(it 

io JK 

<Cel I [ p e \v x \ 2 dxdt + Cel [ [ ( P e fdxdt 

JO JK Jo Jk 

< C(t,K)e* +CS f [ ( P y +1 dxdt ^ £ ^ 0. 

io Jk 

Using (gTgp) and (j4"7?5j) . we obtain 

?$(p e , m e ) t + H^{p\ m t ) x are compact in W^' pl (R+) for some 1< pi < 2. (4.77) 
Furthermore by (f2T2T|) - ([2728]) . lemma HE and (14371) . we nave: 

r]f(p e ,m € ) t + Hf(p e ,m € ) x are uniformly bounded in Lf 3 c (M^_) for p 3 > 2. (4.78) 

where P3 = 7 + 1 > 2 when 7 G (1,3], and P3 = ^j-q > 2 when 7 > 3. By interpolation 
we conclude the proof of the lemma [H 



5 Proof of theorem 11.2 



From lemmas [H we have verified the conditions (i)-(iii) of theorem 11.31 for the sequence 
of solutions (p e ,m e ). Using theorem 11.31 there exists a subsequence (p € ,m e ) and a pair 
of measurable functions (p, m) such that 

{p e ,m e ) ->• (p,m), a.ee-^0. (5.79) 

It is easy to check that (p, m) is a finite-energy entropy solution (p, m) to the Cauchy 
problem (jl.5p with initial data {po,Po u o) for the isentropic Euler equations with 7 > |. 
It achieves the proof of theorem 11.21 
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